Ising spin glass transition in magnetic field out of mean-field 
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The spin-glass transition in external magnetic field is studied both in and out of the limit of 
validity of mean-field theory on a diluted one dimensional chain of Ising spins where exchange 
bonds occur with a probability decaying as the inverse power of the distance. Varying the power in 
this long-range model corresponds, in a one-to-one relationship, to change the dimension in spin- 
glass short-range models. Evidence for a spin-glass transition in magnetic field is found also for 
systems whose equivalent dimension is below the upper critical dimension at zero magnetic field. 
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Introduction - Even though 30 years have passed since 
the spin-glass (SG) phase in presence of an external mag- 
netic field has been characterized in mean-field theory [l| , 
its existence in realistic finite-dimensional systems is not 
yet an established issue. In most common (Heisenberg- 
like) amorphous magnets, e.g., AgMn, CuMn, and AuFe, 
a SG phase has been detected also in presence of an exter- 
nal field [2| ■ In mean-field theory of vectorial spin-glasses 
this transition is expected along the so-called Gabay- 
Toulouse line [3j. In Ising- like materials, instead, like 
Fe x Mni_ x Ti03, it is still a matter of debate whether 
or not a SG phase occurs when the system is embed- 
ded in a magnetic field 0, 0|- Irreversible phenomena 
are, actually, detected in experiments as the temperature 
is lowered: the separation of zero-field cooled and field- 
cooled susceptibilities (or magnetizations) and the rapid 
increase of characteristic relaxation times. In zero field 
these are the signatures of a thermodynamic transition, 
but in some recent AC measurements in a magnetic field 
their magnitude tends to depend sensitively on fre- 
quency and they are interpreted as pertaining to a glassy 
dynamic arrest, rather than to a true thermodynamic 
transition. According to this, the SG features measured 
in a field would be artifacts of being out of equilibrium, 
similarly to what happens in the structural glass transi- 
tion, in which the liquid glass former falls out of equilib- 
rium at some low T when its structural relaxation time 
becomes longer than the observation time and it vitrifies 
into an amorphous solid. 

The Replica Symmetry Breaking (RSB) theory, hold- 
ing in the mean-field regime for spin-glasses, predicts a 
thermodynamic transition in magnetic field h at a finite 
temperature [Bj]. In this framework, a transition line, 
called Almeida-Thouless (AT) Q line, can be identified 
in the T — h plane between a paramagnetic and a spin 
glass phase. At sufficiently low dimensions (i.e. below 
the lower critical dimension, Z?l) the transition disap- 
pears. The value of Dl is not known, but it is quite 
possible that in a field is higher than for h — (as it 
happens for a ferromagnet in a random field) . There are 
some numerical evidences (and analytic results) support- 



ing D\ = 2.5 at zero field. For h > 0, some arguments 
suggested but recently Temesvari [f| argued that 

the AT line can be continued below D = 6. In the droplet 
theory, instead, no transition is predicted to remain as 
soon as an infinitesimal field is switched on, indepen- 
dently from the value of D. A crossover length £d(h,T) 
is introduced Q, beyond which the SG phase is destroyed 
by the field. The predictions of TNT scenario [H should 
be similar to those of the droplet model. Extensive nu- 
merical works on the Edwards- Anderson model in 4D 
and 3D yielded evidence both in favor of a transition in 
field 0, and against it 11, 12, [HI]. Unfortunately, 
finite size corrections are very strong in the presence of 
an external field and it is hard to say whether these sim- 
ulations were really testing the thermodynamic limit. To 
overcome this problem we use a recently introduced SG 
model [3] , which can be simulated very efficiently, and a 
new data analysis, which should be less sensitive to finite 
size effects. We report numerical evidences for a ther- 
modynamic phase transition in the presence of external 
fields also in systems for which the mean-field approxi- 
mation is not correct. 

The model - We investigate a one dimensional chain of 



L Ising spins (er^ = ±1) whose Hamiltonian reads [14 1 
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The quenched random couplings are independent and 
identically distributed random variables taking a non 
zero value with a probability decaying with the distance 



between spins Oi and cr_- 
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Non-zero couplings take value ±1 with equal probability. 
We use periodic boundary conditions and a z = 6 av- 
erage coordination number [231 ]. The random field hi is 
Gaussian distributed with zero average and standard de- 
viation h 24j . We will denote the average over quenched 
disorder, both bonds and fields, by an overline. The uni- 
versality class depends on the value of p. For p > 1 it 
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turns out to be equal to the one of the fully connected 
version of the model 15j , where bonds are Gaussian dis- 
tributed with zero mean and a variance depending on the 
distance as J?- oc r~ p . As p varies, the model displays 
different behaviors |l4(: for p < pu = 4/3, the mean-field 
(MF) approximation is exact, while for p > pu, it breaks 
down because of infrared divergences (IRD). The value 
pu = 4/3 corresponds to the upper critical dimension of 
short-range spin-glasses in absence of an external mag- 
netic field (-Dtj = 6). At p > pl = 2 no finite temperature 
transition occurs, even for h — [16]. A relationship be- 
tween p and the dimension D of short-range models can 
be expressed as p = 1 + 2/D which is exact at D\j = 6 
(pu =4/3) and approximated as D < D v . Indeed, the 
lower critical dimension Z? L — 2.5 JL7| corresponds to 
p ~ 1.8, which is 10% less than /?l- We note that in the 
ferromagnetic (ordered) Ising case on the same kind of 
lattices a simple theoretical argument tells us that the 
value of pl is 2 for h — and 1.5 in a field. 

Simulations details and data analysis - To study the 
critical behavior of the model in external field we sim- 
ulated two replicas a> : using the parallel tempering 
(PT) algorithm Field values are h = 0,0.1,0.2,0.3 
for p = 0,1.2,1.4 and h = 0,0.1,0.15 for p = 1.5. We 
used sizes up to L = 2 14 spins for h — and up to 
L = 2 12 for h > 0. The number of samples is between 
32000 and 64000 for all sizes. Thermalization is guaran- 
teed by the logarithmic binning (in base 2) of data in MC 
steps until at least the last two points coincide. The 
presence of SG long range order can be deduced from the 
study of the four-point correlation function 
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and its Fourier transform C(k) |25j. Indeed, both the SG 
susceptibility 

Xsg = C(0) (4) 
and the so-called second-moment correlation length [l9[ 
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diverge at the critical temperature in the thermodynamic 
limit. For finite (but large enough) systems, the following 
scaling laws hold in the MF regime (1 < p < 4/3) 



£1/3 X V" v " ~ C >J ' L^/ 3 

with v = l/(p— 1), and in the IRD regime (p > 4/3) 

^=x(LHT-T e )), i=i(Lt(T-T c j). (7) 

with 2 — r] = p — 1 . Unfortunately finite size corrections 
to the above scaling laws are known to be very large, 
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FIG. 1: C(k)- 1 vs. x /sin(fc/2)/7r at p = 1.5 (IRD regime), 
h = 0.1, T = 2.1 and L = 2 K , with k = 6, . . . , 12. Upper inset: 
quadratic fit to L — 2 12 data, excluding k — 0. Lower inset: 
comparison between C(0) -1 and its extrapolated valued A. 



especially in the presence of an external field. It is very 
important to understand these finite size effects (FSE) 
and try to keep them under control. 

In the main panel of Fig. [1] we plot 1/C(k) versus 
[sin(fc/2)/7r]''~ 1 for an interesting case (IRD regime with 
field). The choice of the variables is dictated by the fact 
that for L — > oo and T > T c the propagator on the lattice 
at small wave numbers should behave like 



C{k)~ 



A + B[sm{k/2)Y 



(8) 



with Xsg = l/A and £ cx (B/^)V0>-1) = (B X s e ) 1/{p ' 1) ■ 
In other words, A(L = oo, T) goes to zero at T c , while 
B{L = oo,T = T c ) stays finite. 

We observe in Fig. [T] that the largest FSE in C(k) are 
in k = 0, which is the data used for estimating \sg- More- 
over FSE for k > have an opposite sign with respect 
to those in k = (cf. lower inset) and consequently £, 
which is a function of C(0) /C(2ir/L), may be strongly af- 
fected. The reason why FSE become smaller increasing 
k is simple: they are more evident in the large x tail of 
C(x) and, thus, at small k in C(k). Moreover, the large 
x part of C{x) strongly depends on (q), which is known 
to have large sample-to-sample fluctuations in a field and 
FSE due to negative overlaps which should disappear in 
the thermodynamic limit. 

With the aim of reducing FSE, we introduce a method 
for estimating T c using C(k) data with k > 0. We 
fit C(k)^ 1 by a quadratic function A + By + Cy 2 with 
y = [sin(fc/2)/7r]' 3 ~ 1 : the goodness of such a fit can be ap- 
preciated in the upper inset of Fig.[TJ As long as T > T c , 



we expect limx^oo A(L, T) 



> 0: the lower inset 



in Fig. Q] shows size dependence of C(0) _1 and A(L,T), 
having compatible thermodynamic limits. 

In the main panel of Fig. [5] we show the best fitting 
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FIG. 2: Plot of A(L, T) vs. T at p = 1.5, h = 0.1. Sizes are 
L = 2 K , with « = 7, . . . , 12. Inset: T C (L) vs. L" ' 28 . 



FIG. 4: Plot of Xs g /L 0A vs. T at p = 1.4, /i = 0.1. Sizes are 
L = 2 K , with k = 6, . . . , 12. Inset: £/£ vs. T. 
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FIG. 3: Plot of Xsg/i°' 5 vs. T at p = 1.5, /i = 0.1. Sizes are 
L = 2 K , with k = 6, . . . , 12. Inset: £/L vs. T. 



FIG. 5: Plot of Xsg/^ 1/3 vs. T at p = 1.2, h = 0.2. Sizes are 
L = 2 K , with ac = 6, . . . , 12. Inset: £/L" /3 vs. T. 



parameter A(L,T) for p = 1.5 and /i = 0.1. For each 
size we compute the temperature T C (L) by solving the 
equation A(L, T C (L)) = (in this way only A > data 
are used, which are the most reliable). Finally, we es- 
timate T c = limi-xx) T C (L) (inset of Fig. [2|) and obtain 
T c = 1.46(3). The T C (L) scaling in L~ x l v has an expo- 
nent —0.28, in good agreement with 1/v = 0.25(3) for 
the h — case [141 ]. On the same data (p — 1.5, h = 0.1) 
the analysis of the crossing points of Xsg/L 2 ~ v and £/L, 
cf. Eq. j7]), is shown in Fig. [3l yielding no evidence for 
a phase transition. The most natural explanation is the 
presence of corrections to scaling laws Eq. (0 . 

The case p = 1.4 and h — 0.1 provides a still more use- 
ful comparison. Our method of analysis returns a critical 
temperature T c — 1.67(7). Fig. [4] shows Xsg/L 2 ~ v and 
£/L vs. T: crossings are present, but the curves seem to 
merge for T < 1.5 and a precise determination of T c is 



practically unfeasible. For p — 1.2, h = .2, cf. Fig.[5j the 
estimate based on the scaling of x sg ~ L 1 / 3 - Eq.© - 
yields T c — 1.67(3), while the ^/L"/ 3 curves do not show 
any crossing for T > 1.2. Since the transition is known 
to be there in MF, this behavior of £ is clearly caused by 
large FSE. 

Numerical values for the estimates of T c obtained with 
the two methods are reported in Table [J and look com- 
patible. It is clear that for large p our method works 
better. As p is decreased, this new estimate becomes 
poorer, because the scaling exponent p — 1 [cf. Eq.©] is 
too small to yield a robust extrapolation of A(L,T). 

Discussion of experimental results - A possible objec- 
tion to the presence of the SG transition (supported by 
our results) is that in experiments on Ising-likc SG no AT 
line was detected. Here wc consider, in particular, the 
most recent experiments on Feo.55Mno.4sTi03 4], where 
the AC susceptibilities were very accurately measured in 
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TABLE I: Estimates of T c : column 5 from Eqs.lJBHB and 
column 6 from the extrapolation of A(L,T) by Eq.([8]). 
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transition line in the diluted power-law decaying interac- 
tion Ising SG chain at all values of the power analyzed, 
including values corresponding to short-range SG models 
below the upper critical dimension. The behaviour be- 
low the AT line may change with the dimension. We are 
presently studying this possibility. 

These AT lines were not found in the study of the fully 
connected version performed in Ref . [2lJ , nor in Ref . [22j ] 
where a similar diluted model was simulated [2f|. There, 
T c was estimated by using the scaling properties of £,/L. 
As we have shown, this quantity suffers of strong FSE. 
We put forward an alternative method to discriminate 
between a pure paramagnetic phase at all temperatures 
and a finite temperature spin-glass transition. One of the 
advantages of this method is that it mainly uses data at 
T > T c , leading to more reliable results. 

For what concerns three dimensional real systems, we 
hint that the magnitude of the external fields used in ex- 
periments up to now might be too large to firmly rule out 
the presence of an AT transition line. We suggest a range 
of fields (h < 1000 Oe) where the transition should take 
place in Feo.ssMno^TiOa and we hope this may stimu- 
late further experimental investigations. 

This work has been partially supported by MEC, con- 
tracts FIS2006-08533-C03 and FIS2007-60977. Part of 
simulations were performed in the BIFI cluster. 



FIG. 6: Relative decrease of the critical temperature with 
increasing field (p = 1.5, h = 0,0.1,0.15,0.2) versus relative 
decrease of x*> the ZFC susceptibility at T c (h = 0). The 
dashed curve is a guide for the eye. 

the presence of an external magnetic field. In order to 
relate external fields in our model to those used in ex- 
periments we look how much the zero-field-cooled (ZFC) 
susceptibility at temperature T c (h = 0), x*> decreases 
as h is increased. In Fig. [6] we plot T c (h)/T c (0) versus 
X*(^)/x*(0) m our model for p = 1.5. In experiments 
on Feo.55Mno.45Ti03[Ij with external fields of magnitude 
100 Oe, 300 Oe, 1000 Oe and 3000 Oe one has, respec- 
tively, x*O)/x*(0) = 0.98, 0.94, 0.84, 0.75. These ra- 
tios, cf. Fig. [6l suggest that a SG transition is unlikely 
to be experimentally observed above h — 1000 Oe. 

Increasing p, and/or considering J„ ^ 0, the critical 
field decreases. The p — 1.5 model considered above 
is approximately equivalent to a short-range system in 
D = 4. Therefore, in order to detect, or rule out, a SG 
phase in D = 3, it becomes even more important to work 
at small fields. The observation that the fields used in 
experiments on Ising-like materials are maybe too large 
to see a SG phase is in agreement with the results of 
Petit et al., who studied both Ising-like and Heisenberg- 
like spin glass samples 0. 

Conclusions - In conclusion, by using a new method 
of data analysis, we have been able to identify an AT 
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